Conjugate 18cm OH Satellite Lines at a Cosmological Distance 
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We have detected the two 18cm OH satellite lines from the z ~ 0.247 source PKS 1413+135, the 1720 MHz 
line in emission and the 1612 MHz line in absorption. The 1720 MHz luminosity is -Loh ~ 354Lq, more than 
an order of magnitude larger than that of any other known 1720 MHz maser. The profiles of the two satellite lines 
are conjugate, implying that they arise in the same gas. This allows us to test for any changes in the values of 
fundamental constants, without being affected by systematic uncertainties arising from relative motions between 
the gas clouds in which the different lines arise. Our data constrain changes in G = g P [a. 2 /y] 1 ' 849 , where 
y = m e /m p ; we find AG/G — 2.2 ± 3.8 x 10 _J , consistent with no changes in a, g P and y. 

PACS numbers: 98.80.Es,06.20Jr,33.20.Bx,98.58.-w 



INTRODUCTION 

In recent times, much interest has centred on the possibility 
that fundamental "constants" such as the fine structure con- 
stant a might vary with cosmic time (e.g. |l|, |2], |3[). Many 
theoretical models, including Kaluza-Klein theories and su- 
perstring theories, predict spatio-temporal variation of these 
constants. However, terrestrial experiments have so far shown 
no evidence for such changes. For example, the strongest 
presently available constraints on changes in a arise from iso- 
topic abundances measured in the Oklo natural fission reactor, 
which give the constraint Aa/a < 1.2 x 10~ 7 on fractional 
changes in the fine structure constant |4]. Uzan |5] provides 
an excellent review of current experimental and observational 
constraints on changes in a and other fundamental constants. 

While terrestrial studies such as the Oklo experiment have a 
high sensitivity, these measurements only probe a small frac- 
tion of the age of the Universe. For example, the Oklo reactor 
operated about 1.8 Gyr ago, less than a sixth of the age of the 
Universe 13111 . Thus, current terrestrial experiments cannot 
rule out earlier changes in the fundamental constants, quite 
possible if these changes are non-monotonic in nature. 

Astrophysical studies of redshifted spectral lines provide a 
powerful probe of putative changes in fundamental constants 
over large fractions of the age of the Universe (e.g. |lLy,|g,|2l|). 
Perhaps the most interesting of these estimates arise from 
the new "many-multiplet" method (jj, |8J, which has been 
used in conjunction with Keck telescope spectra to obtain 
Aa/a = (-0.54 ± 0.12) x 10~ 5 over the redshift range 
0.2 < z < 3.7 |9] (see, however, 1 10]). On the other hand, 
the many-multiplet method was recently applied to Very Large 
Telescope (VLT) data to obtain Aa/a = (-0.6 ±0.6) x 10~ 6 
over the redshift range 0.4 < z < 2.3 |3], in direct con- 
flict with the Keck results. Further, many theoretical analy- 
ses predict that changes in a should be accompanied by much 



larger changes in other constants, such as the ratio of elec- 
tron mass to proton mass m e /m p (e.g. II lL il'M ). How- 
ever, fractional changes in m e /m p have been found to be 
less than (—0.5 ± 3.6) x 10 -5 over a similar redshift range 

(0 < z < 2.75) H. 

Given the importance of changes in fundamental constants 
for theoretical physics, it is important that this be tested by 
independent techniques. This is especially true since the 
same technique (i.e. the many-multiplet method) has yielded 
strongly divergent results when applied to data from differ- 
ent telescopes, suggesting that systematic errors play a signif- 
icant role (e.g. 1 14]). We have recently (|7], hereafter Paper 
I; see also I150 developed a new technique to simultaneously 
constrain changes in a, the ratio of electron mass to proton 
mass y = m e /m p and the proton gyromagnetic ratio g p , using 
the multiple 18cm OH lines. This makes use of the fact that 
these lines arise from two very different physical mechanisms, 
Lambda-doubling and hyperfine splitting 1 16], and thus have 
very different dependences on a, y and g p . If g p is assumed 
to remain unchanged (e.g. |2, 6]), observations of all four OH 
18cm transitions in a single cosmologically distant galaxy can 
be used to simultaneously estimate any changes in y and a. 
Further, since these lines arise from the same species, a com- 
parison of the OH column density obtained from each line can 
be used to test whether the lines arise from the same gas. This 
is often an issue in other techniques, since lines from different 
species are used in the comparison, giving rise to systematic 
uncertainties arising from the possibility of relative motions 
between the gas clouds in which the different lines arise. 

Only the "main" 18cm OH lines (at frequencies of ~ 1667 
and 1665 MHz) have hitherto been detected at cosmologi- 
cal distances OHO; the satellite lines (at ~ 1612 and 
1720 MHz), which are critical to the above test, have not been 
found until now. We report here the first detection of the satel- 
lite OH lines at cosmologically significant distances, from the 



0.247 source PKS1413+135. 



SPECTRA AND RESULTS 

Spectra of the redshifted OH 1612 and 1720 MHz lines 
from PKS1413+135 were obtained with the Westerbork Syn- 
thesis Radio Telescope (WSRT) in June 2003, while the 1665 
and 1667 MHz transitions were observed with the Giant Me- 
trewave Radio Telescope (GMRT) in October 2001. The top 
three panels of Fig. 1 show the 1720, 1612 and 1667 MHz 
spectra, respectively, while the bottom panel shows the sum 
of the 1612 and 1720 MHz optical depth profiles; the velocity 
resolution is ~ 1.1 km/s for the satellite lines and ^3.5 km/s 
for the 1667MHz line. Note that PKS1413+135 is unresolved 
by the GMRT and WSRT beams; the optical depths plotted 
here are the ratio of line depth to the total flux density for the 
different transitions. The RMS noise in each of the satellite 
spectra is ~ 0.002 (per 1.1 km/s channel) while that in the 
1667 MHz spectrum is ~ 0.0008 (per 3.5 km/s channel), in 
units of optical depth. 

The most striking feature of Fig. 1 is that the 1612 and 
1720 MHz lines are conjugate to each other; this can be clearly 
seen in Fig.^D], where the sum of the 1612 and 1720 MHz 
optical depth profiles is consistent with noise. We note that 
the weak feature seen at ~ —20 km/s in Fig.[nD] is not sta- 
tistically significant (< 3er, even after smoothing the profile 
to increase the signal-to-noise ratio). The peaks of the two 
profiles also arise at the same redshift, within the error bars. 

Such conjugate behavior of the satellite lines has been seen 
earlier in extra-galactic sources 1 20, 21] and essentially arises 
due to competition between two decay routes to the 2n 3 / 2 (J = 
3/2) ground state, after the molecules have been pumped from 
the ground to the higher excited rotational states 12211 . The 
pumping can take place either due to collisions or far-infrared 
(FIR) radiation. If the last step of the cascade is the 119/i 
transition 2n 3/2 (J = 5/2) -> 2n 3/2 (J = 3/2), the 1720 MHz 
line is inverted and the 1612 MHz line anti-inverted; con- 
versely, if the final step is the 79/j decay 2^/2 (J = 1/2) 
— + 2n 3 / 2 (J = 3/2), it results in 1612 MHz inversion and 
1720 MHz anti-inversion iQ. The fact that the 1720 MHz 
line is seen in emission and the 1612 MHz line in absorp- 
tion in PKS1413+135 implies that the former route dominates 
and gives the constraint Nqh/AV < 10 15 cm -2 (km/s) -1 
||20j], where iVoH is the OH column density and AV, the 
velocity width; in our case, AV ~ 20 km/s, implying that 
Woh J$ 2 x 10 16 cm~ 2 . Further, all pumping mecha- 
nisms require that the 119/i transition be optically thick for 
the 1720 MHz line to be inverted; this corresponds to the 
lower limit Nrm/AV > 1 x 10 14 cm" 2 (km/s)" 1 £3l, i.e. 



iVoH > 2 x 10 cm in the present case. We thus ob- 
tain the range 2 x 10 15 < iV OH < 2 x 10 16 cm" 2 for the 
OH column density. At lower -/Voh, the 1 19/i line is optically 
thin and the 1720 MHz line is not inverted; at higher AToh, 
the inversion switches to the 1612 MHz line. Implications for 
other physical conditions in the absorbing/emitting cloud will 



be discussed elsewhere. 

The 1720 MHz profile implies an OH luminosity Loh ~ 
354L Q , well into the OH megamaser range 12311 . This is by far 
the brightest known 1720 MHz maser line, more than an order 
of magnitude more luminous than Arp 220 {Lou ~ 12£©), 
the only other known 1720 MHz megamaser 1 24]. Interest- 
ingly enough, PKS1413+135 is an excellent candidate for 
main line OH megamaser emission, on the basis of its high 
FIR flux I23J I25B : in fact, the empirical relation, log Loh = 
(1.2±0.1)logL F iR-(H-7±1.2) (in solar units) Hj predicts 
a main line megamaser luminosity of Loh ~ 300L Q , similar 
to that obtained in the 1720 MHz line. It is unclear whether 
this is merely a coincidence, especially given the spread in the 
relation between Lqh and Lfir 12511 . 



CONSTRAINING THE VARIATION OF FUNDAMENTAL 
CONSTANTS 

Crucially, the conjugate behavior of the satellite lines guar- 
antees that they arise from the same gas. However, the main 
line profiles are quite different, suggesting that this absorp- 
tion occurs at a different spatial location. The main lines (the 
1665 MHz line is not shown here) peak at z = 0.24671, co- 
incident in velocity with other molecular species (such as CO 
and HCO+) that have been detected in PKS1413+135 Efill . 
They also show a "tail" at positive velocities, with a total ve- 
locity spread of ~ 15 km/s. On the other hand, while the two 
satellite lines have a similar velocity width (~ 20 km/s), they 
show blueshifted, edge-brightened profiles, with the absorp- 
tion/emission peaking ~ —15 km/s away from z — 0.24671. 
Further, while the satellite features are seen to extend to v = 0, 
the main lines are not detected at the velocity corresponding 
to the peak of the satellite lines. 

Most of the astrophysical techniques used to estimate 
changes in fundamental constants (e.g. |3J |6j |9J) involve a 
comparison between the redshifts of spectral lines of different 
species, i.e. the assumption that these different species have 
no velocity offsets between them. In fact, even a compari- 
son between lines of the same species (e.g. 1 13, 27]) does not 
rule out the above systematic uncertainty, as the different lines 
might well be excited in different regions of the gas cloud. In 
the present case, the very different shapes of the profiles of 
the "main" and "satellite" line profiles implies that one cannot 
compare the peak redshifts of all four lines to simultaneously 
constrain changes in y = m e /m p and a. However, the conju- 
gate behavior of the OH satellite lines implies that they arise 
from precisely the same gas; we hence have a unique situation 
where systematic velocity offsets between two spectral lines 
are not an issue. Again, the satellite line frequencies have 
very different dependences on the "constants" a,y = m e /m p 
and the proton g-factor g p ; this can be clearly seen from equa- 
tions (9) and (11) in Paper 1 (note that the sum of satellite 
line frequencies is equal to the sum of main line frequencies 
in the OH radical). A comparison between the measured 1612 
and 1720 MHz redshifts thus allows one to measure changes 
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FIG. 1: OH optical depth r v/s heliocentric velocity (relative to z — 0.24671). [A] WSRT 1.1 km/s resolution 1720 MHz spectrum. [B] 
WSRT 1.1 km/s resolution 1612 MHz spectrum. [C] GMRT 3.5 km/s resolution 1667 MHz spectrum. [D] Sum of WSRT 1612 and 1720 MHz 
spectra. The spectrum is consistent with noise. 



in the quantity G = g p [or/y] . In particular, if z s and Zd 
are the redshifts derived from the sum and difference of the 
measured 1612 and 1720 MHz line frequencies, respectively, 
it can be shown that AG/G = [z s — Zd]/{1 + z), where z is 
the average of z s and Zd- 

As mentioned earlier, the peaks of the two satellite profiles 
arise at the same redshift, within our observing resolution of 
~ 4.8 kHz. However, since the two profiles are conjugate, 
we can improve the accuracy of the comparison between the 
peak redshifts by fitting a template profile to each spectrum. 
Note that the same template is fit to both spectra and this pro- 
cedure is carried out solely for the purpose of comparing the 
peak redshifts. We use a double-Gaussian template for this 
purpose, to account for both the narrow peak as well as the 
extended shoulder. Several of these two-component fits were 
found to be statistically indistinguishable; in essence, the pro- 
file cannot be uniquely decomposed into two Gaussians as the 



signal-to-noise ratio of the shoulder is quite low. Using dif- 
ferent templates allowed us to estimate the systematic error in 
the measured redshifts, arising due to the choice of template. 
In all cases, the redshifts of the peak were found to agree to 
within 2.8 kHz, for both the 1720 MHz and 1612 MHz lines. 
We hence use this figure as our estimate of the error in the red- 
shifts (including both systematic and random errors). We note 
that attempts were also made to fit a single Gaussian compo- 
nent to the narrow peak in both satellite line profiles; again, 
the peak redshifts were found to be in agreement within an 
error of 2.8 kHz. 

The sum of the frequencies of the peaks of the satellite pro- 
files is then measured to be 2673. 356±0. 004 MHz, while their 
difference is 86.8733 ± 0.004 MHz. In both cases, the errors 
have been added in quadrature - it should be emphasized that 
this is a conservative estimate in the case of the difference of 
the frequencies, since the systematic error in the choice of the 



template cancels out here, to first order. The above frequen- 
cies correspond to redshifts z a = 0.246658 ± 0.0000015 and 
z<± = 0.24663 ± 0.000046, respectively, in agreement within 
the error bars. This gives the limit AG/G = 2.2 ± 3.8 x 10~ 5 
on changes in the quantity G from z ~ 0.247 to the present 
epoch. Again, we emphasize that the errors on AG/G are 
dominated by the error on the difference in frequencies and the 
latter is certainly lower than 2.8 kHz, since the systematic er- 
ror from the fitting procedure cancels to first order; we choose, 
however, to remain conservative in our error estimates. 

We note that the strong dependence of the quantity G on the 
fine structure constant a (G oc a 3698 ) implies that the obser- 
vations are extremely sensitive to changes in a. If we assume, 
as is often done, (e.g. |6]) that neither g p nor y varies with cos- 
mological epoch, the above constraint on AG/G translates to 
Aa/a = 0.6 ± 1.0 x 10~ 5 . Similarly, G also has a strong de- 
pendence on the electron-proton mass ratio y = m e /m p (al- 
though weaker than that on a), with G oc y -1 - 849 ; we obtain 
the constraint Ay/y = (—1.2 ± 2.0) x 10~ 5 , assuming that a 
and g p remain unchanged. This is similar to the best current 
limits on changes in y from observations of the Lyman and 
Werner band lines of molecular hydrogen 1 13], although the 
latter constraints, of course, probe a larger redshift range. We 
also note that the changes in a and y are of opposite sign and it 
would require something of a conspiracy for both these quan- 
tities to vary, while leaving G unchanged. The data towards 
PKS1413+135 are thus consistent with the different constants 
remaining constant from z ~ 0.247 to today. 

The relative sensitivity of the present technique adds to its 
importance. While results from the many-multiplet method 
(e.g. 1 3, 9]) on fractional changes in a have lower errors at the 
present time, this method requires the use of many absorption 
systems, both to cancel systematic effects and to increase sen- 
sitivity. Further, as mentioned earlier, the conflicting results 
from analyses based on Keck and VLT spectra suggest that 
unknown systematics dominate the errors here. On the other 
hand, our present limits are derived from a single system and 
provide an entirely independent constraint on changes in a, y 
and g p . It should also be noted that the present satellite spec- 
tra have a resolution of only ~ 1.1 km/s; higher resolution 
spectroscopy in both lines will allow far tighter constraints on 
changes in G = g p [a 2 /y] 1 ' 849 . For example, observations of 
PKS1413+135 with the Square Kilometer Array O, a next 
generation radio telescope, will be able to detect fractional 
changes AG/G ~ few x 10~ 7 . This implies a sensitivity of 
~ 10~ 7 to changes in Aa/a, comparable to those obtained 
with the Oklo reactor experiment. 

Finally, the 1720 MHz inversion arises due to an over- 
population of the F = 2 level of the ground state relative 
to the F = 1 level, as transitions from the F = 3 level 
of the J = 5/2 state to the latter level are forbidden. This 
should also give rise to an over-population of the above F = 3 
level, relative to the J = 5/2, F = 2 level (by ~ 50 %; 
12911 ') and could result in the inversion (and anti-inversion) of 
the corresponding J = 5/2 satellite lines (at rest frequen- 
cies of — 6049 MHz and 6016 MHz). A detection of these 



lines would allow one to independently constrain changes in 
all three constants a, y and g p 13011 : these observations are 
presently being carried out. 

In summary, we report the first detection of the 18cm OH 
satellite lines at cosmological distances. The 1720 MHz line 
is found in emission, with a luminosity larger by more than an 
order of magnitude than any other known 1720 MHz maser; 
the 1612 MHz line is in absorption, as are the other ground 
state OH lines. The satellite line profiles are found to be con- 
jugate, in that the sum of their optical depths is consistent with 
noise; this implies that the two lines arise in precisely the same 
gas. A comparison between the satellite line redshifts then 
yields constraints on G = g p [a 2 /y] 1849 , a combination of 
the fine structure constant a, the electron-proton mass ratio 
y = m e /m p and the proton gyromagnetic ratio g p . We find 
AG/G = 2.2 ± 3.8 x 10~ 5 , consistent with the different con- 
stants remaining unchanged from z ~ 0.247 (a lookback time 
of ~ 2.9 Gyr) to the present epoch. If y and g p are assumed to 
remain constant, we obtain Aa/a = 0.6 ± 1.0 x 10~ 5 , over 
the redshift range < z < 0.247. 
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